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Abstract. We present a detailed microscopic derivation for a non-Mar kovian master 
equation for a driven two-state system interacting with a general structured reservoir. 
The master equation is derived using the timc-convolutionless projection operator 
technique in the limit of weak coupling between the two-state quantum system and its 
environment. We briefly discuss the Markov approximation, the secular approximation 
and their validity. 
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1. Introduction 

All quantum systems are open, i.e., they interact with their environment. This leads 
to dissipation and to a loss of quantum coherence due to an irreversible flow of energy 
and/or information from the system to its environment [HE]. In some cases, however, 
it is possible that the properties of the environment modify the dynamics of dissipation 
and decoherence. In particular, reservoir memory effects can damp or even momentarily 
reverse the direction of the flow of energy and information, and the quantum system 
can regain some information that was previously lost to the environment [H H] . 
The properties of the environment are described by the spectral density function 
J{uj). When the spectral density is structured, i.e., J{ui) varies considerably with 
frequency, environmental memory effects can become significant. Physical examples 
of environments with a structured spectral density include photonic band-gap materials 
and leaky optical cavities [5J [6] . 

When the reservoir memory effects considerably modify the dynamics of the quantum 
system, the dynamical semigroup property is violated and a non-Markovian description 
of the dynamics is called for [7]. A systematic method for deriving a non-Markovian 
master equation is given by projection operator techniques [8J [9]. In general, a non- 
Markovian master equation contains an integration over the whole past history of the 
system of interest. However, in the limit of weak coupling between the system and the 
environment one may use the time-convolutionless (TCL) projection operator technique 
to obtain a simpler local-in-time non-Markovian master equation [TO] . 
In this paper we derive the non-Markovian master equation for a damped driven two- 
state system in the limit of weak coupling between the quantum system of interest and 
its environment. We assume the environment to be described by a generic structured 
spectral density. A microscopic master equation for the damped driven two state system 
in the Markovian regime has been derived by Tanas and Kowalewska-Kudlaszyk for the 
special case of a Lorentzian environment [11]. The non-Markovian dynamics of such a 
system has been studied by Budini, using a phenomeno logical master equation [12] . We 
present a microscopic derivation for a non-Markovian master equation. The microscopic 
master equation is useful when we wish to understand the microscopic processes that 
determine the dynamics of the quantum system. 

2. Derivation of the non-Markovian master equation 

2.1. The model 

We consider a two-level atom with an energy separation uja {h~ = 1) driven by a 
monochromatic laser of frequency ujl almost resonant with the frequency of the atom 
\b>A ~ Wl\ to a- The two-level atom interacts with a quantised electromagnetic field 
modelled by an infinite chain of quantum harmonic oscillators of frequencies {ujk}. For 
convenience we study the dynamics of such system in a frame rotating at the laser 
frequency ujl- In the rotating frame the free Hamiltonians of the system of interest and 
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of the environment, and the interaction Hamiltonian are 

H s = ^(Aa z + na x ), (1) 
H E = J2 UJ ka k a{, (2) 

k 

Ei = ^2g k [exp{ko L t)a + a k + exp(-iu; L t)cr_4], (3) 

k 

respectively. Here (T x ,y,z are the usual Pauli operators and a± are the raising and lowering 
operators. The operators are given in the atomic basis, i.e., an orthonormal basis 
{|e), \g)} determined by the excited and ground state of the two-level atom, fl is the 
Rabi frequency of the laser and g k is a coupling term describing the interaction between 
the two-level atom and the k-th quantum harmonic oscillator. We also defined the 
detuning A = oja — ojl- 

The system Hamiltonian ([1]) and the interaction Hamiltonian ([3]) are both given in the 
rotating wave approximation, i.e., we have neglected non-energy conserving terms. 



2.2. The TCL2 master equation 

According to the time-convolutionless projection operator technique a general master 
equation to second order in coupling between the system of interest and its environment 
is 

d ,. x /•* 



dt p{t) = -J Q ds tr E [#/(*), [H T (s), p{t) ® p E ]\ , (4) 



where 



Hj(t) = ex V [it(H s + HE^Hjexpi-itiHs + H E )\ (5) 



is the interaction Hamiltonian in the interaction picture. We have assumed factorized 
initial conditions for the total density matrix ptot{fy — p(0)®p_b(0). Furthermore, in the 
limit of weak coupling, the state of the environment remains stationary and therefore 
PE(t) = pe(0) = pE for all times t. 



2.3. Transformation to the eigenbasis 

For convenience we go to the eigenbasis of the free system Hamiltonian H g, i.e., 
orthonormal basis defined by 

|V>±) = ±^= (VI ± sin 6\e) + V1T sin %)) , (6) 

where 9 = arctan(A/fi). In the eigenbasis representation the system Hamiltonian is 
diagonal, H s = (u/2)a z , where uj = \J A 2 + uj 1 is the energy bias between the eigenstates 
and a bar over an operator indicates that the operator is given in the eigenbasis. The 
interaction Hamiltonian in the interaction picture becomes 

Hi(t) = £ Si(t)®Ei(t), (7) 

i=l,2 
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where 



Si(t) = 6XP ^ L ^ [exp(-io;*)(A - u)a_ + exp(iut)(A + uo)a + + Vta z ] ,(8) 



2u 



k 

and S 2 (t) = S\(t) and E 2 (t) = E\(t). 



(9) 



2.4- The environment 

We assume that the environment is in a thermal state p E = exp (— H^/hsT), where 
Z is the partition function and k% is the Boltzmann constant. We will consider the case 
when temperature T = and only a single reservoir correlation function is non-zero: 

tT E [E i {t 2 )E j {ti)p E ]=0, i^2,j^l, 

tT E [E 2 (t 2 )E 1 (t 1 )p E ] = 5>xp[k^(t 2 - hM. (10) 



In the limit of a continuum of modes of the environment 
J2gt -> / dajJ(u), 



where we introduce the generic spectral density J (to). Combining equations (j4]), ([7]), 



poj) and pj we find 

^p(t) = Ijdr | du/ J^Oe-^^Ct-r)^)^^) 

- S , i(t)S , 2 (t-r)p(t)]+ J ff.c, (12) 
where H.c. is an abbreviation for Hermitian conjugation. At this point we define 

r e (t) = jf* dr y dcj' J(w') exp[i(c; L - w' + fw)t], £ G {-1, 0, 1}. (13) 

Recall that the detuning is small, i.e., |A| <C u>a,^l- Furthermore, in physically 
realistic situations the Rabi frequency is also small, i.e. Q <C ua,^l- Therefore 
uj = VA 2 + fi 2 C and T±i(t) ~ To(^) = r(t). In the following we decompose 
the quantity T(t) into real and imaginary parts, i.e. T(t) = (l/2)j(t) + i\(t). 
Keeping this in mind we obtain form equation (112p the master equation for the damped 
driven two-state system in the Schrodinger picture 



dt 



p{t)= -i[H s + H L ]+C 2 +1 (t) 



a-p(t)a+ - ]-{a + a-,p(t)} 



+ chit) 



a+p(t)a- - -{(j_a+,p(t)} 



v z p{t)o- z - ~{a z a z ,p{t)} 



+ cut) 

- C_C 7(£) [a + p{t)a z + a z p(t)a_ 
+ C + C j(t) [o-p(t)o z + a z p(t)a + 
+ Cl 1 {t){a + p{t)a + + cr_p{t)dJ\ 
'lit) 



+ C 



-{a x ,p(t)} + i\(t)[a xi p(t)} 



(14) 
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where the Lamb shift Hamiltonian is 

H L = \{t)[C\a_a + + C 2 + a_a_ + C^a% (15) 

and 

c± = c ° = aJ- (16) 

2.5. Markov approximation 

The reservoir correlation time Tc is defined by the properties of the environment. If 
Tc is small compared to the relaxation time of the system tr, the reservoir memory 
effects are negligible and one can make the Markov approximation. When the Markov 
approximation is valid one can extend the upper limit of integration in equation ([TBI 
from t to infinity. In the Markovian limit the decay rates j(t) become time independent 
positive constants 7 = lim^ 7(t). 

We note that in the limit of weak approximation between the system of interest and 
its environment the first non-Markovian corrections do not change the operatorial form 
of the Markovian master equation. Instead the first non-Markovian corrections only 
affect the decay rates by making them time- independent. For times t = 0{tc) the time- 
dependent decay rates deviate from their stationary Markovian values and it is during 
this timescale that non-Markovian reservoir memory effects take place. 



2.6. The secular approximation 

The typical timescale for the damped driven two-state system is t$ = uj~ x . When 
the typical timescale is much smaller than the relaxation time r# one can make the 
secular approximation (rotating wave approximation) on the master equation. The 
secular approximation comprises of neglecting terms that, in the interaction picture, 
oscillate rapidly with frequency r^ 1 . In the case considered in this paper the secular 
approximation amounts to neglecting the last four lines of the master equation (1141) . 
When the secular approximation is valid the master equation consists of a coherent part 
describing unitary evolution generated by Hs + Hl, and of a part with three terms in 
a Lindblad-type form, i.e., in a Lindblad form with time- dependent coefficients. The 
Lindblad-type terms describe quantum jumps that disrupt the unitary evolution of the 
system. 

When 7(i) > one can describe the quantum jumps in the framework of Monte Carlo 
wave functions [13]. In this case the quantum jumps lead to irreversible dissipation and 
decoherence of the quantum system. If, however, j(t) < for some period of time, a 
more general description, namely the non-Markovian quantum jump (NMQJ) method 
is called for [3]. For time interval when the decay rate is negative the direction of 
the quantum jumps is reversed. Reversed quantum jumps reverse dissipative processes 
and allow the quantum system to regain some of the information and/or energy that 
was previously lost to the environment. The system can recohere and reconstruct 
superposition states. The NMQJ description explains how memory effects affects the 
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dynamics of a quantum system described by a local-in-time non-Markovian master 
equation. 

3. Conclusions 

We have derived a local-in-time non-Markovian master equation for a damped driven 
two-state system interacting with a general structured environment using the time- 
convolutionless projection operator technique to second order in coupling between the 
system and its environment. We assumed factorized initial conditions and a stationary 
reservoir. To keep master equation (fl4j) as general as possible we did not perform the 
secular approximation. 

We found that to second order in perturbation theory the operatorial form of the non- 
Markovian master equation is the same as the operatorial form of the corresponding 
Markovian master equation. The important difference that arises from the first non- 
Markovian corrections is that, for timescales specified by the reservoir memory time, 
the decay rates become time-dependent. 

The specific form of the decay rates can be computed once the environment specific 
spectral density function has been specified. Master equation (HM lays the basis for a 
systematic study of non-Markovian dynamics of the damped driven two-state system 
embedded in different types of reservoirs. 
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